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Excited state spectra are calculated using lattice QCD for baryons that can be formed from 
u, d and s quarks, namely the A'', A, A, E, H and SI families of baryons. Baryonic operators 
are constructed from continuum operators that transform as irreducible representations of SU(3)f 
symmetry for flavor, SU(4) symmetry for Dirac spins of quarks and 0(3) symmetry for orbital 
angular momenta. Covariant derivatives are used to realize orbital angular momenta. Using the 
operators, we calculate matrices of correlation functions in order to extract excited states. The 
resulting lattice spectra have bands of baryonic states with well-defined total spins up to J = |. 
Each state can be assigned a dominant flavor symmetry and the counting of states of each flavor and 
spin reflects SU{6) x 0(3) symmetry for the lowest negative-parity and positive-parity bands. States 
with strong hybrid content are identified through the dominance of chromo-magnetic operators. 



I. INTRODUCTION 

The spectra of baryon resonances have been a focus of 
study experimentally and theoretically, both for particles 
composed of the light u and d quarks, and those contain- 
ing one or more of the heavy c and b quarks. There is 
increasing activity aimed at understanding the spectra of 
particles containing one or more s quarks, so-called hy- 
peron physics, e.g., [ll. At present the knowledge of the 
S and families is particularly limited with only a few 
states experimentally established, and scant knowledge 
as to their properties [2] . Two recent lattice studies have 
focused on the lowest few states of [3] and A j4j . 

A quantitative description of the spectra of baryons 
that can be constructed from the u, d and s quarks 
is important for a number of reasons. Firstly, to ex- 
plore how the effective degrees of freedom that describe 
the hadron spectra change as the masses of the quarks 
are changed. Secondly, to advance our understanding of 
QCD in regimes where the hyperons play a crucial role, 
such as in the physics of the early universe and core- 
collapse in supernovae. Finally, to address the long-term 
goal of extracting baryon spectra from lattice QCD. 

In this paper, we present lattice QCD calculations of 
the excited-state spectra of baryons that can be con- 
structed from u, d and s quarks: the A'^, A, A, S, S 
and fl families of baryons. Our calculations exploit a 
menu of methods that we have developed and are key 
to our studies of excited-state spectra: the use of an 
anisotropic, clover action for the generation of gauge con- 
figurations [3 |6] , the construction of operators that re- 
spect the symmetries of the lattice yet which retain a 
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memory of their continuum analogues [7| , the use of "dis- 
tillation" [5] to efficiently compute the correlation func- 
tions between those operators, and finally the application 
of the variational method [51 [TD] , exploiting the eigenvec- 
tors of matrices of correlation functions to determine the 
spins and flavors of the extracted energies. In two earlier 
works, Refs. [Hill], we applied these methods to compute 
the nucleon and A spectra, with the spins of the states 
identified clearly; these works revealed spectra at least 
as rich as the quark model, with suggestions of "hybrid" 
states in which the gluons played an important structural 
role. Here we expand upon these works by considering 
baryons containing not only the light u and d quarks, but 
also one or more s quarks. 

The layout of the remainder of the paper is as follows. 
In the next section, we describe details of the lattices 
used, and outline the construction of the lattice interpo- 



lating operators. In Sec. II B we recall our procedure for 
analyzing the hadron correlation functions, and for iden- 
tifying the continuum spins of the states. Our results 



are presented in Sec. Ill and a summary of the work is 
presented in Sec. |IV[ 



II. 



COMPUTATIONAL METHODS AND 
DETAILS 
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A significant challenge in determining the excited-state 
spectra is to obtain a sufficient number of energy lev- 
els in order to extract their patterns of energies, spins 
and flavors before statistical noise overcomes the sig- 
nal. This requires accurately resolving the behavior of 
hadron correlation functions at short temporal separa- 
tions. A computationally efficient calculation is obtained 
through the use of an anisotropic action, with a finer 
temporal lattice spacing than that used in the spatial 
directions, enabling correlation functions to be resolved 
over several time slices while preserving a sufficient spa- 
tial volume. The lattice action used in this work, as well 
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with maximum accessible values of orbital angular mo- 
mentum, L, of 0,1, and 2, respectively. Tables [lT| and 
nil show the patterns of quantum numbers of operators 
based on non-relativistic spins. Not listed are operator 
constructions based on the lower components of Dirac 
spinors (relativistic quark spins), although they outnum- 
ber the non-relativistic operators in the set used. 

Three non-relativistic quark spins are mixed- 
symmetric for S = ^ and symmetric for 5* ~ |. 
Orbital angular momenta are mixed-symmetric for 
L = 1 and symmetric, mixed-symmetric or antisymmet- 
ric for L = 2. The combination of flavor, spin, and space 
symmetries must be symmetric to yield an operator that 
is antisymmetric when color is included, in accord with 
the Pauli Principle. We list in the tables the spin-parity, 
, of the distinct combinations that are allowed. The 
classification of operators according to SU{3)f flavor 
symmetry will be exploited as a means of identifying the 
dominant flavor structure of the states in the extracted 
spectrum. 

The operators constructed from non-relativistic quark 
spinors have the spin, parity and flavor quantum num- 
bers allowed by SU{6) x 0(3) symmetry. These quan- 
tum numbers occur in definite patterns that are indicated 
by the bold numbers in Tables [ll] and [ITU there are 
iVi(J), Ns{J) and Nio{J) operators for SU(3)f singlet, 
octet and decuplet symmetries, respectively, with total 
angular momentum J. If the states in the spectra were 
to correspond to broken SU{6) x 0(3) symmetry, then 
there would be patterns with the same numbers of states 
and the same quantum numbers. We will compare our 
spectra with the SU{6) x 0(3) patterns given in Tables [n] 
andUm 

For d = 2, the mixed-symmetry combination D]^Li 
corresponds to the commutator of two covariant deriva- 
tives, producing a chromo-magnetic field, as explained in 
Ref. [11] . We call these hybrid operators because of their 
essential gluonic content. The hybrid operators based on 
non-relativistic quark spinors are listed in Table [IV| The 
patterns of states expected to be created by such hybrid 
operators correspond to the numbers Mi (J), Ms{J) and 
Afio(J). 

The final step in the construction of the interpolating 
operators is the subduction to the irreducible representa- 
tions (irreps) of the cubic group. See Ref. [Tj for details. 
The number of operators used in each lattice irrep is the 
same for both positive and negative parities, which are 
denoted by the subscripts g and u, respectively, and the 
operators are classified according to the flavor irreps of 
SU{3)f, as shown in Table |v| 

In order to reduce the complexity of the analysis, 
some relativistic operators are omitted. We include all 
non-relativistic operators based on the upper-component 
Dirac spinors. For example, for E, 68 operators in irrep 
H are used, with 37 having flavor octet symmetry and 31 
having flavor decuplet symmetry. This set omits 22 rela- 
tivistic operators but is sufficient to determine the spec- 
trum. Because we omit some relativistic operators, and 



-0.0743 -0.0743 702 1.00 0.3593(7) 500 7 56 
-0.0808 -0.0743 524 1.15 0.3200(7) 500 7 56 
-0.0840 -0.0743 391 1.39 0.2951(22) 479 8 56 

TABLE I: Parameters of the 16'' x 128 lattices and propagators 
used in this work. The pion mass in MeV, and the number of 
configurations are listed, as well as the number of time-sources 
and the number of distillation vectors A^vccs. 

as the method used to tune the parameters of the ac- 
tion, can be found in Refs. [3 [6]. To summarize, we 
use improved gauge and fermion actions with two mass- 
degenerate light quarks of mass mi and a strange quark 
of mass TTis. We employ 16^ x 128 lattices having spa- 
tial lattice spacing ~ 0.123 fm, and a renormalized 
anisotropy, the ratio of the spatial and temporal lattice 
spacings, of ^ « 3.5. The calculation is performed at 
three values of the light-quark masses, corresponding to 
pion masses of 391, 524 and 702 MeV. The 702 MeV pion 
mass corresponds to the SU{3)p flavor-symmetric point. 
Some details of the lattices are summarized in Table [Tj 

A. Baryon Interpolating operators 

The construction of the baryon interpolating operators 
is described in Ref. [7| . Our starting point is to construct 
a set of continuum baryon interpolating operators, which 
we express symbolically as, 

where the factors F, S and D describe the flavor, Dirac 
spin and derivative structure of the operators, respec- 
tively, and the subscripts, S^, denote the corresponding 
symmetry representations with respect to permutations. 
The Dirac spin factor, (5'^^)^^, represents the combina- 
tion of three Dirac-spinor quark fields with overall an- 
gular momentum S and parity Ps, with permutational 
symmetry Eg, while n labels different constructions. The 
spatial factor, -D^'^u' expressed in terms of gauge- 
covariant derivatives acting on the three quark fields, 
where [d] denotes the number of derivatives, S^j is the 
permutation symmetry, and the derivatives are combined 
so as to transform as angular momentum L. Finally, the 
factors are combined to yield an interpolator of definite 
spin and parity, which we label . 

The construction of operators is guided by symmetry 
considerations. Based upon four-component Dirac-spinor 
quark fields combined with SU(3) fiavor plus angular 
momentum, the symmetry of the full set of operators is 
SU{12)xO{3). An important subset of operators is based 
on non-relativistic quark spins, by which we mean opera- 
tors based on the upper two components of Dirac-spinor 
quark fields. In this study, we realize angular momenta 
by including up to two covariant derivatives, d = 0, 1, 2, 
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8f 1 1 

2 1 

2 ^ 


1 - 3 - 

2 2 

1 - 3 - 5 - 

2 2 2 


iV8(J) 


2 2 1 


lOp 1 1 


1 - 3 - 

2 2 


iVio(J) 


110 


If \ 1 


1 - 3 - 

2 2 


iVi(J) 


110 



TABLE II: Allowed spin-parity patterns for one-derivative op- 
erators based on non-relativistic quark spinors {SU{Q) x 0(3) 
symmetry). For each flavor, the quark spin S and orbital 
angular momentum L are listed followed by the allowed 
values. The total number of operators is listed as for flavor 
octets, A''io for flavor decuplets and A''i for flavor singlets. 



we include the "hybrid" operators introduced in Ref. [TT] , 
the number of operators differs from Ref. [7 . 

The spectra for the baryons made from u and d quarks, 
i.e., N and A, have been explored in earlier works using 
the same ensembles and operator constructions we will 
employ here, firstly in a calculation of the spectra for 
both parities but with the "hybrid" operators excluded 
from the basis [7], and then of the positive-parity spectra 
using all operators, including the "hybrid" operators [TT]. 
In this work we consider excited states of both parities 
for baryons that can be formed from u, d and s quarks. 

B. Correlator Analysis 

The variational method we use in our analysis involves 
the computation of a matrix of correlation functions, 

Ct,{t^tf-U) = {0\O,{tf)O]{U)\Q), (2) 

for operators i and j that lie in a given cubic irrep, A. 
With the large operator basis introduced above, it is 
essential to have an efficient computational method for 
computing the baryon correlation functions. We use dis- 
tillation [8^, which provides a smearing function whilst 
enabling the full matrix of correlators to be constructed 
for all the operators at both source and sink. Following 
our earlier work using these lattices, we employ iV = 56 
eigenvectors of the gauge-covariant Laplacian when con- 
structing the distillation operator, and compute correla- 
tion functions from N^^^^^ time sources, as listed in Ta- 
ble H 

Our fitting strategy follows exactly that outlined in 
Refs. and jllj . In summary, we solve the generalized 
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TABLE III: Allowed spin-parity patterns for two-derivative 
operators based on non-relativistic quark spinors {SU{6) ® 
0(3) symmetry). For 8f, two distinct operators with different 
internal symmetries are allowed for L=0 combined with S=| 
and for L=2 combined with S=|. Hybrid operators are listed 
separately in Table|lV] The total number of operators for each 
spin is listed as A'^s for flavor octets, TVio for flavor decuplets 
and A''i for flavor singlets. 

eigenvalue equation, 

C{t)v„{t, to) = X„{t, to)C{to)v„{t, to), (3) 

and thereby obtain the masses from the principal corre- 
lators, Xn{t, to), n = 1 . . . dim(C), by fitting to, 

Xnit,to) - (1 - A„)e-™"(*-*") + A„e-""(*"*«). (4) 

The corresponding eigenvectors, Vn{t,to), are aligned 
with those at reference time slice, t = tj-^i, and they pro- 
vide information as to the optimal interpolating operator 
for the state n, namely J^i'^n^l- Figure [l] shows fits of 
principal correlators for the ground state and eight ex- 
cited states in the S spectrum for irrep Hg, where each 

state is identified as having spin-parity = f^- These 
states will be discussed later when we show how their 
spins and flavors can be identified and when we show 
that they form part of a pattern of states within the E 
spectrum. 

As in Ref. , we make extensive use of the operator 
"overlap factors", Z" = (n|O||0) that occur in the spec- 



4 



14 - 



Hg 'Hi,,, HI:- 
State 

I .... I ■ . 



10 eo 
t 




1.0 - 



X V^Joi^ i 4/ 1 5: ' m = 6,4'077(3Z 



Hg 

- State 2 



_L 



_L 



10 so 

t 




Hg 
- State 5 



_L 



h i Hg 

i state 7 

: ■ ■ ■ I I 




" h ^/Nd»(- ^ 6/ 1 5; ■ m - 6.4377(28 



1.0 - 



Hg 

state 13 



'■'•••lilllil 



10 



10 



1.0 - 




FIG. 1: Fits to principal correlators for nine states in irrep Hg that are identified as J = Fits are obtained using Eq. (j4j|. 
For plotting, we divide out the first exponential factor, thus, the plots show values of e™"''"'"' An(t) as the data points. Lines 
show the fits according to the form e'""(*~*oU„(t) = l~An + A„ e"'™""™-")'*"*"), with to = 9. The fits approach the constant 
value, f — j4i,, for large t. The bands indicate the fits with the central fit values along with one-standard-deviation uncertainties. 



tral decomposition of the matrices of correlation func- 
tions, 



(5) 



Using the orthogonality condition for the eigenvectors 



v^^C{tQ)v^ — (5n,mi the overlap factors are related to the 
eigenvectors through Zf = yST^e'""*''/^ y^*Cj,{to). 

The baryon operators used in this work are constructed 
in two stages. The first stage is to construct operators 
in the continuum that have definite spin quantum num- 
bers. In principle, these operators would produce a ma- 
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TABLE IV: Allowed spin-parity patterns for hybrid two- 
derivative operators based on non-relativistic quark spinors, 
following Ref. [11]. The operators correspond to the combi- 
nation of three quarks in a color octet with a gluonic field, 
G, to make a color singlet as [{qqq)saGsa] i ■ The gluon field 
has spin-parity 1^. The patterns of states expected to be cre- 
ated through the coupling of a chromo-magnetic gluon field 
coupled to quark fields in a color octet correspond to the num- 
bers AIs for flavor octets, Mio for flavor decuplets and Mi for 
flavor singlets. The counting of operators is the same as in 
Table In] but the operators have reversed parity. 
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TABLE V: For each 5L''(3)f symmetry, isospin / and 
strangeness S, the numbers of operators are shown in the last 
three columns for each lattice irrep Gi, H and G2. The con- 
structions use operators with up to two derivatives. The same 
numbers of operators are used for both positive and negative 
parities. 



trix of correlation functions that is orthogonal, i.e., pro- 
portional to 5jji, where J and J' are the spins of the 
source and sink operators. In the second stage of the 
construction, the operators are subduced to the irreps 
of the cubic group so that they can be used on a finite 
lattice. We observe that after the subduction, there re- 
mains a remarkable degree of rotational symmetry in the 
matrices of correlation functions that are obtained for 
baryons. They exhibit approximately the orthogonality 
property that holds in the continuum. An example of 
this is shown in Fig. [2j where the S correlator matrix. 
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FIG. 2: Magnitudes of the elements of the correlator matrix, 
Cij I sjCiiCjj, at time-slice 5 are shown for E at m,r = 391 
MeV according to the darkness scale at the bottom. The 
blocks along the diagonal correspond to operators with the 
following spin-flavor combinations: |(10f), |(8f), |(10f), 
|(8f), |(10f) and |(8f). 



Cij{t), is shown at t=5 for irrep Hg, after normalizing 
each operator so that the diagonal elements are equal to 
one. The matrix indices range over 68 operators as fol- 



lows: 19 J 
J 



|, lOp operators; 22 J : 



|, IOf operators; 12 J = |, 8f operators; 3 J 



|, 8f operators; 9 

_ 7 

IOf operator; and 3 J = |, 8f operators. The matrix is 
approximately block-diagonal with respect to J. Similar 
block-diagonal matrices of correlation functions are ob- 
served within other irreps and they stem from the con- 
struction of operators starting from definite continuum 
spins. These approximately orthogonal matrices are key 
to the identification of the continuum spin of states. 

Although flavor symmetry remains a broken symmetry 
in the continuum because the strange quark in our calcu- 
lations is heavier than the u and d quarks, the correlator 
matrix is approximately block diagonal with respect to 
flavor. As shown in Ref. p^, which uses the same lat- 
tices, the breaking of S't/(3)F symmetry is weak also for 
the mesons. When the operators of a given flavor sym- 
metry are dominant, we exploit that feature to identify 
the dominant flavor composition of a state, and when 
hybrid operators play a substantial role we use that to 
identify states with strong hybrid content. Figure [3] illus- 
trates how such identifications are made for states in the 
Hg irrep of the E. The plot shows the overlaps, Z", of 
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operators that create each state, with operators grouped 
by their continuum spin and flavor as indicated in the 
column labels at the top. The columns also are labeled 
by whether the operators are based on relativistic (R) or 
non-relativistic (NR) Dirac spinors, or hybrid construc- 
tions (h) or non- hybrid ones (nh). The operator overlaps 
provide a reasonable identification of the spin, flavor and 
hybrid content of the S states in irrep Hg of the cubic 
group. In particular, the identifications of states with 
spin-parity | , and flavor irreps Sp or IOf, are noted in 
the caption. Similar identifications are made to extract 
the patterns of other states in our spectra. 



III. RESULTS 

In the following, we present the spectra in units of the 
f2 mass determined on the corresponding ensemble. The 
spectra of TV, A, A, E, S and Q, families of baryons for 
spins up to I and both parities are shown at two different 
pion masses in Figures |4] and [Sj The dominant flavor 
irrep is indicated by color: blue for Sp, yellow for IOf, 
and beige for If. At the S'[/(3)F-symmetric point, the 
spectra are shown in Fig. |6j now classified according to 
their SU{'ii)F flavor irrep: Sp, lOp or Ip. Results for the 
positive parity N and A spectra were presented earlier 
in ref. [11]. 

A general feature of the lattice spectra is that there 
are bands of excited states with alternating parities for 
each family of baryons: the lowest band of states for each 
parity is shown inside the slanted boxes. In the lowest 
negative-parity band, the numbers of states of each spin 
and flavor agree with the expectations shown in Table [III 
the table that shows the quantum numbers allowed by 
SU{Q) X 0(3) symmetry for one-derivative operators. 

In the lowest positive-parity band, the numbers of 
states for each spin and flavor agree with the expectations 
shown in Table |III[ the table that shows the quantum 
numbers allowed by SU{(!>) x 0(3) symmetry for two- 
derivative operators. While at the SU{?>)f symmetric 
point the counting involves a single flavor irrep, for the 
A, E and S, the correct numbers of states are obtained 
from summing the numbers for each of the two flavor ir- 
reps involved. This agreement between spectra and the 
expectations based on non-relativistic quark spins pro- 
vides a clear signature of SU{6) x 0(3) symmetry in the 
spectra. 

As noted in the caption of Fig. [3) the overlaps of 
the different operators classified according to their flavor 
structure provides a means of identifying the dominant 
flavor composition of the states in the spectra. Although 
the mixing of flavors increases as one moves away from 
the S'[/(3)i?-symmetric point, the dominant flavor iden- 
tifications are shown by the different colors used in the 
plots of spectra. It is clear from our analysis that for 
the mixed-flavor states, the A, E and 5, we find spec- 
tra that exhibit the multiplicities expected from exact 
SU{Z)f flavor symmetry for each of the flavor-identified 



multiplets. In summary, the baryon spectra are remark- 
ably consistent with the SU{&) x 0(3) expectations. It 
is this symmetry that is the basis for the quark model 
and the lowest bands of lattice states include the same 
quantum numbers that occur in the quark model |13H15j . 

In addition to the three-quark states that correspond 
to the quark model, a number of states are identified by 
their strong overlaps with hybrid interpolating operators, 
implying that they have a strong hybrid content. These 
states, in which the gluons play a substantive role, are 
shown for positive parity by symbols with thick borders 
in Figs.jijjsjandle] The two J = E states in Fig. [i] 
with strong hybrid content and masses near l.Gmn are 
Hg states 16 and 17 in Fig. [Sl Note that the states with 
strong hybrid content generally are at high energy, typi- 
cally about O.Ttoq, or more, above the ground state. 

In this work we have no three-derivative operators. Be- 
cause of that, and because the relativistic operators gen- 
erally play an important role at higher energies, it is not 
meaningful to interpret the observed states at higher en- 
ergies in terms of SU{6) 0(3) symmetry. However, one 
sees qualitative similarity of the higher bands in Figs. |4] 
and [5] to the corresponding bands found at the flavor- 
symmetric point in Fig. |6] 

It also is not meaningful to identify states with strong 
hybrid content in our negative-parity spectra and none 
are shown. Signals for states with significant hybrid con- 
tent exist in negative-parity states at high energy, but 
they are based on the relativistic operators, whereas the 
three-derivative, non-relativistic operators that are ab- 
sent may be equally, or more, important. Without a clear 
understanding of the relative importance of all the rele- 
vant operators, we cannot identify negative-parity states 
with strong hybrid content. 

The baryon spectrum does not admit the "spin-parity 
exotics" that provide a useful indication of hybrid states 
in the meson spectrum, and indeed we observe that hy- 
brid operators can have a significant, albeit not domi- 
nant, contribution to many states. Some examples can 
be found in Fig. [sj notably the lowest- lying spin- 1 decu- 
plet. 

The patterns and multiplicities of positive-parity states 
with strong hybrid content can be compared with the 
expectations based on non-relativistic quark spins that 



are listed in Table IV At the flavor symmetric point with 

m^r — 702 MeV, which is shown in Fig. [G] and at 

— 524 MeV, which is shown in Fig. [s) one sees all the 

positive-parity hybrid states corresponding to Table |IV| 

At TOtt — 391 MeV, which is shown in Fig. |4l most of 

1 + 

the hybrid states are seen. One decuplet, J — ^ hybrid 
state is missing for E. An octet, >/ = f ^ state and a 

decuplet, = f state are missing for ^. These baryons 
involve mixings of 8p and lOp flavor symmetries and 
their spectra are particularly dense when subduced to the 
lattice irreps. For example, an Hg state must be found 
for each >/ = |, | and | state in the spectrum. The S, 
J = I state with strong hybrid content was found as the 
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FIG. 3: "Matrix" plot of values of operator overlaps, , for state n and operator i, normalized according to 



so that 



the largest overlap across all states for a given operator is unity. For each of the n = to 21 states of E in the Hg irrep, the 
magnitude of each operator's overlap is shown for = 391 MeV. Darker pixels indicate larger values of the operator overlaps 
as in Fig. [5] Column labels indicate non-relativistic (NR) and relativistic (R) operators, as well as hybrid (h) and non-hybrid 
(nh) operators. In addition the flavor irrep is indicated as (IOf) for decuplet or (8f) for octet and continuum spins of the 
operators are shown by |, | and |. State 0, the ground state, and excited states 2, 8 and 15 are identified as = | states 
with dominant decuplet flavor symmetry. States 1, 4, 5, 7 and 13 are identified as — excited states with dominant octet 
symmetry. States 16 and 17 are — f ^ excited states with strong hybrid content and dominant octet flavor symmetry. 



29 state in the Hg spectrum and was the highest state 
determined. The fact that a few of the states observed at 
higher values of the pion mass are not found at our lowest 
value of m^r may be because not all states have been 
determined. Overall, the lowest states with strong hybrid 
content are in reasonable accord with the expectations 
based on Table |IVl 

The excited states in the lowest bands of negative- 
parity are particularly well determined: Fig. [7] shows the 
ones that are created predominantly by flavor-octet oper- 
ators and Fig. [8] shows the ones that are created predom- 
inantly by flavor decuplet and singlet operators. Note 
that we show the data for the mass of state n, namely 
TTin, in physical units that are obtained from the formula, 
TOn = 1672.45 where mo latt is the il mass on the 

ensemble. Thus, the physical ft mass is used to set the 
scale. The patterns of these states are very similar for 
the different baryons. For the flavor-octet states, there 
are two ^ states, two | states and one | state. For 
each baryon, the energy increases with spin J and the 



highest energy is about 300 MeV above the lowest en- 
ergy, independent of the baryon. For the flavor decuplet 
case, there is one 1 state and one | state with about 
70 to 100 MeV splitting, Sio being an exception. The 
flavor singlet case has the same pattern except that the 
energies are lower and the splitting is larger. 



We note that there are several important limitations of 
the present study. They have been discussed in Ref. [7] 
and we conclude with a brief summary of them. The 
16^x128 lattice used is small, with spatial dimensions 
of about 1.9 fm on a side. The pion masses used are 
significantly larger than the physical mass. No operators 
that efficiently couple onto scattering states (e.g., ttN) 
are included. Studies of the resonances that correspond 
to the three-quark states will require improvements that 
overcome each of these limitations. 
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FIG. 4: Results for baryon excited states using the ensemble with = 391 MeV are shown versus J^. Colors are used to 
display the flavor symmetry of dominant operators as follows: blue for 8f; beige for If; yellow for lOp. Symbols with thick 
border lines indicate states with strong hybrid content. Calculations are for a 16^ x 128 lattice. The lowest bands of positive- 
and negative-parity states are highlighted within slanted boxes. The eight excited states of E, with ^ that are shown 

within a slanted box, are Hg states 1, 2, 4, 5, 7, 8, 13 and 15. Fits for the same states are shown in Fig. [irand identifications 
of their spins and flavors are noted in Fig. [3] 




FIG. 5: Results for baryon excited states using the ensemble with = 524 MeV are shown versus . Symbols are as 
described in Fig. |4] 
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FIG. 6: Results for baryon excited states in flavor irreps 8f, IOf and If obtained using the flavor-symmetric point, with m,r 
= 702 MeV, are shown versus . Symbols are as described in Fig. IZ] 
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FIG. 7: The lowest negative-parity states that are flavor- 
octet are shown for = 391 MeV. 



FIG. 8: The lowest negative-parity states that are flavor- 
singlet (beige) and decuplet (yellow) are shown for m^r = 391 
MeV. 



IV. SUMMARY 

This work presents results for baryons based on lattice 
QCD using the 16^x128 anisotropic lattices that were 
developed in Ref. [5] . Excited state spectra are calculated 
for baryons that can be formed from u, d and s quarks, 
namely the N, A, A, S, S and $7 families of baryons, for 
two pion masses, 391 MeV, 524 MeV, and at the SU{i) p- 
symmetric point corresponding to a pion mass of 702 
MeV. 

The interpolating operators used incorporate covariant 
derivatives in combinations that correspond to angular- 
momentum quantum numbers L = 0, 1 and 2. The an- 
gular momenta are combined with quark spins to build 
operators that transform according to good total angu- 
lar momentum, J, in the continuum. As noted in earlier 
works, approximate rotational symmetry is realized at 
the scale of hadrons, enabling us to identify reliably the 
spins in the spectrum up to J 
a single lattice spacing. 

The operators we have employed are classified accord- 
ing to the irreducible representations of SU(3)f flavor. 
At the pion masses used, the SU(3)i? symmetry is broken 
only weakly and states in the spectra can be identified 
as being created predominantly by operators of definite 
flavor symmetry 8i?, lOp or If. 

We find bands of states with alternating parities and 
increasing energies. Each state has a well-defined spin 
and generally a dominant flavor content can be identi- 
fled. The number of non-hybrid states of each spin and 
flavor in the lowest-energy bands is in agreement with the 



I from calculations at 



expectations based on weakly broken SU{6) ^0{3) sym- 
metry. These states correspond to the quantum numbers 
of the quark model. 

Chromo-magnetic operators are used to identify states 
that have strong hybrid content. Usually these states 
are at higher masses, about 0.7mn, or more, above the 
lowest non-hybrid states. There is reasonable agreement 
between the number of positive-parity states with strong 
hybrid content and the expectations of Table [TV] that are 
based on non-relativistic quark spins, although a few of 
the expected states are not found at the lowest pion mass. 

With the inclusion into our basis of multi-hadron oper- 
ators, which couple efficiently onto multi-hadron scatter- 
ing states, we expect to find an increased number of levels 
in the spectrum. As demonstrated in Ref. |16| , using the 
technique of moving frames where the total momentum of 
the system is nonzero, the increased number of levels al- 
lows for the extensive mapping of the energy dependence 
of scattering amplitudes, and hence, the determination of 
resonances. The prospect of determining the properties 
of resonances provides a strong motivation for continued 
work on the spectra of baryons. 
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